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ABSTRACT 

Observations of the starburst galaxy, M82, have been made with the VLA in its A- 
configuration at 15 GHz and MERLIN at 5 GHz enabling a spectral analysis of the 
compact radio structure on a scale of < 0.1" (1.6 pc). Crucial to these observations 
was the inclusion of the Pie Town VLBA antenna, which increased the resolution of 
the VLA observations by a factor of ^2. A number of the weaker sources are shown 
to have thermal spectra and are identified as Hil regions with emission measures 
~10 7 cm" 6 pc. Some of the sources appear to be optically thick at 5 GHz implying 
even higher emission measures of ^10 8 cm -6 pc. The number of compact radio sources 
in M82 whose origin has been determined is now 46, of which 30 are supernova related 
and the remaining 16 are Hn regions. An additional 15 sources are noted, but have 
yet to be identified, meaning that the total number of compact sources in M82 is at 
least 61. Also, it is shown that the distribution of Hn regions is correlated with the 
large-scale ionised gas distribution, but is different from the distribution of supernova 
remnants. In addition, the brightest Hn region at (B1950) 09 h 51 m 42.21 s +69° 54' 
59.2" shows a spectral index gradient across its resolved structure which we attribute 
to the source becoming optically thick towards its centre. 
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1 INTRODUCTION 

The radio emission from star-forming galaxies can be ei- 
ther of a thermal (e.g. free-free) or non-thermal (e.g. syn- 
chrotron radiation) origin. The non-thermal radiation orig- 
inates from relativistic electrons which have been acceler- 
ated by supernovae resulting from the deaths of massive 
stars and the thermal emission is from free-free emission 
from gas which has been ionised by hot, young stars. There- 
fore, it is obviously of interest to investigate the relative 
contributions of these processes to the emission which we 
observe from nearby star-forming galaxies. A separation of 
the non-thermal and thermal components of radio emission 
from the archet ypal starburst galaxy, M8 2, has already been 
carried out by Allen & Kronberg (1999) on angular scales 
of 2" and they identified several complexes of ionised gas. 
However, 2" is not sufficient resolution to separate the most 
compact radio sources from the more diffuse radio struc- 
ture. In this paper we compare observations of M82 made 
with the VLA in its A-configuration at 15 GHz, incorpo- 



rating the VLBA Pie Town antenna, with MERLIN 5 GHz 
observations. The maps made from these observations have 
a resolution < 0.1" which extends the analysis of the most 
compact radio structure into the thermal regime. Thus, we 
are now able to comment on the relative contribution of 
thermal and non-thermal processes to the radio emission on 
scales of order 1 pc. For consistency with previous publica- 
tions, we assume a distance of 3.2 Mpc to M82 throughout 
the paper (Burbidge, Burbidge & Rubin 1964). 



Spectral information on 26 of the co mpact radio sources 
M82 has already be en published (Wills et al. 1997; 



Allen & Kronberg 199c). These analyses were prompted 



by 408 MHz MERLIN observations and the radio contin- 
uum spectra of the detected sources were mainly consistent 
with non-thermal synchrotron emission, occasionally with 
a low-frequency turnover which was attributed to free-free 
absorption by ionised gas. However, free-free absorption is 
likely to be responsible for the non-detection of ma ny of 
the remainin g sources at low frequencies. In addition, Wills 



ct al. (1997) (hereafter W97) identified two sources which 
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were inconsistent with a free- free absorbed synchrotron spec- 
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trum, the AGN candi date, 44. 01+59. rjj] and a p ossible Hn 
region at 40.62+56.0. |Allen fc Kronberg (199g| ) (hereafter 
AK98) came to similar conclusions regarding the majority 
of the sources, but additionally identified a flat-spectrum 
(q ~ -0.1, where 5* oc v a ) source at 42.21+59.2, which they 
identified as an Hn region. However, over 50 compact radio 
sources have been identified in the nuclear region of M82, 
therefore until now there has only been spectral informa- 
tion published on around half of the total number, simi- 
l ar to the situation in another nearby starburst, NGC 2146 
( Larchi et al. 2000 ). Also, since previous investigations con- 
centrated on the brightest sources at the longer radio wave- 
lengths, there would have been a bias away from the identifi- 
cation of thermal sources such as Hn regions. This is because 
their electron temperatures are unlikely to exceed 10 4 K and 
hence even if the ionised gas is optically thick, the bright- 
ness temperatures cannot exceed this limit (n.b. 10 4 K = 
0.48 mjy beam -1 with MERLIN at 5 GHz, 50 mas res- 
olution). Therefore, in order to make robust comparisons 
between the numbers of compact radio sources of different 
types in starburst galaxies we must first make more confi- 
dent identifications of a greater fraction of the sources. Fur- 
thermore, with the improved sensitivity of radio observa- 
tions, it is clear that further populations of compact radio 
sources will be discovered in more distant galaxies. Hence 
the need to fully investigate the nature of these sources in 
nearby galaxies. 



2 THE OBSERVATIONS AND DATA 
ANALYSIS 

Observations of M82 were performed with the Very Large 
Array (VLA) of the National Radio Astronomy Observatory, 
including the Pie Town Very Long Baseline Array (VLBA) 
telescope, the data from which were correlated in real-time 
via a fibre-optic link to the VLA site, on the 2nd Decem- 
ber 2000. These observations represent the first stage of a 
programme to image the 15 GHz radio emission from M82 
on a wide range of scales and the combination with B- and 
C-configuation data will be presented in a later paper, along 
with a discussion of the structures and sizes of the compact 
sources. However, the A-configuration data alone represent 
measurements of the radio emission from M82 on scales com- 
parable to those measured using MERLIN at 5 GHz. There- 
fore, we also present observations made with the MERLIN 
interferometer in February 1999. The details of both obser- 
vations are summarised in Table |l[ A number of the param- 
eters listed in Table |l| were chosen to optimise the field of 
view of the observations since the compact sources in M82 
cover a region roughly 50" x 10" in size. Hence, the VLA 
correlator was used in spectral line mode, limiting the usable 
bandwidth to 43.75 MHz, but producing 7 spectral channels 
with widths of 6.25 MHz. This mode results in smearing 
at the 5% level at a distance of 95" from the phase cen- 
tre of the observations (Taylor, Carilli & Perley (1999). We 
deemed this level of smearing to be acceptable for the pur- 
poses of our observations and Figure [l] shows the map of the 



1 All positions within M82 are quoted relative to (B1950) 09 /l 
51 m +69° 54' 



field of view of the observations. Table |l| shows that the two 
sets of observations cover a similar overall range of spatial 
frequencies, although the VLA has more shorter baselines 
and hence greater sensitivity to extended structure. In addi- 
tion, MERLIN has some slightly longer baselines (in terms 
of wavelengths) than the VLA and Pie Town combination. 
Therefore, in order to provide a valid comparison, maps were 
made by limiting both data-sets to a range of 130 kA < L < 
3660 kA, where L = (u 2 + v 2 ) 1 ^ 2 and u and v are the co- 
ordinates of the visibilities in the u-v plane. The lower limit 
to this range means that any structures larger than ~1.6" 
are suppressed. For the purposes of this paper, the word 
'compact' may describe any source which has significant ra- 
dio emission on the scales measured by our data, i.e. less 
than ~1.6" in size. In addition, the density of sampling in 
the u-v plane is substantially different between the VLA and 
MERLIN, since MERLIN is a more sparsely distributed ar- 
ray. Therefore, the relative weighting functions used on the 
two data-sets needed to be different. The final imaging pa- 
rameters used in the comparison are also presented in Ta- 
ble §. 

The MERLIN and VLA data were reduced separately 
and self-calibration was used to remove residual complex 
gain errors after phase-referencing. Therefore, the images 
from each data-set were registered relative to one another 
by 'fixing' the position of the brightest and most compact 
radio sourc e 41.95+57.5, at (B 1950) 09 h 51 m 41.95" +69° 
54' 57.50" (Pedlar et al. 1999). A spectral index map and 
spectral index error map was then produced for the entire 
central region of M82 using the AIPS task COMB. Spectral in- 
dices were only calculated for pixels with brightnesses above 
the 3-cr level in each map as stated in Table 0, using the spec- 
tral index convention stated previously. Spectral index error 
maps were also calculated using the noise levels stated in 
Table |l[ but are not presented. Sufficient signal was present 
in the maps at both frequencies such that spectral index 
maps could be derived for 34 sources. The MERLIN and 
VLA images of these 34 sources are shown in Figure ^ and 
the spectral index maps for 5 of these sources are shown in 
Figure []. For around 75% of the 34 sources in Figure [| there 
exists a good correspondance between the positions of the 
compact sources (within ~10 mas) at 5 GHz and 15 GHz. 
In the remaining ~25% of cases more complex structure is 
present which does not agree as well between frequencies. 
Some possi ble reasons for these discrepancies are discussed 
in Section 3.2. In addition to the 34 sources detected in 



both sets of observations, we detect an additional 5 compact 
sources in the VLA map which have no counterparts in the 
MERLIN map on the same scale. These sources were cho- 
sen by examining the spatial frequency limited VLA 15 GHz 
map and searching for compact sources with emission above 
the 5-cr level. The VLA maps of these sources are shown in 
Figure ^| with the corresponding limits on their two-point 
spectral indices. For reference purposes, Figure jj] contains 
boxes indicating which sources are illustrated in Figure || 
and circles indicating the sources in Figure ^ It should be 
noted that Figure [I] shows a significant amount of 15 GHz 
emission from regions outside of the regions imaged in Fig- 
ures |2] and ^. Most of this emission is suppressed once the 
limitation on the spatial frequency coverage of the 15 GHz 
observations is made, although emission is detected at a ~3- 
a level in these regions. 
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Table 1. Observational setup for the VLA A-configuration (incorporating the VLB A Pie Town antenna) and MERLIN observations of 
M82. 



Array 



VLA (A-configuration) + Pie Town 



MERLIN 



Date 

Observing Frequency 

Bandwidth 
Number of Frequency Channels 
Integration Time per Visibility 
Total Integration Time 
Longest Baseline (kA) 
Shortest Baseline (kA) 
Naturally Weighted Sensitivity 
(/ijy beam -1 ) 
Gaussian u-v taper FWHM (MA) 
CLEAN beam size (mas) 



2nd December 2000 
14.9649 GHz 

43.75 MHz 
7 x 6.25 MHz 
3.33 sees 
12 hours 
3662 
20.80 
50 

4.0 

87x77 & 58° 



5th / 22nd February 1999 
4.546 / 4.866 / 5.186 GHz 
(Multi-Frequency Synthesis) 
15 MHz 
15 X 1 MHz 
4 sees 
24 hours 
3768 
132.1 
35 

2.2 

83x74 @ 42° 



3 THE NATURE OF THE COMPACT 
SOURCES 

By comparing the maps of the compact radio emission in 
M82 at 5 GHz and 15 GHz we have detected 34 sources 
that display compact radio emission at both frequencies as 
discussed in Section and identified an additional 5 com- 
pact sources at 15 GHz, but not at 5 GHz. Of these 39 
sources, 19 have had spectra published by either W97, AK98 
or both. We agree with the original identification of 18 of 
these 19 sources as being of supernova origin and the re- 
maining source as being an Hn region, although the nature 
of the non-thermal sources at 41.95+59.5 and 44.01+59.6 is 



still questionable (e.g. Wills et al. 1999 



McDonald et al 



2001). Therefore, there were 20 sources which remained to 
be reliably identified. In order to do this we examined the 
radio continuum spectrum of each source from 1.4 GHz to 
15 GHz. 



3.1 Spectra of the compact sources 

Since two-point spectral indices only offer a limited amount 
of information regarding the nature of the source it is ob- 
viously of interest to examine the full spectra of as many 
sources as possible. Therefore, we have used lower resolu- 
tion (~0.2" typically) maps at 1.42 GHz (L-band), 4.86 GHz 
(C-band), 8.4 GHz (X-band) and 15 GHz (U-band) and ex- 
amined the flux densities at the positions given by the VLA 
and MERLIN comparison maps. For the L-band measure- 
ments we used the com bined VLA and MERLIN map of 
Pedlar & Muxlow (1995) and for the C-band measurements 



we used an existing VLA A and B-configuration combined 
image (observed on 5th July 1995 and 31st October 1995 in 
A and B-configurations respect ively). For the X-ba nd mea- 



surements we used the data in Huang et al. (1994) for the 
sources available and for the U-band map we used the new 
VLA data and utilised the full spatial frequency range with 
an appropriate weighting for a comparison on a scale of 
~0.2". Since these low-resolution maps contain a substan- 
tial amount of the more diffuse background radio emission 
in M82, an estimate of the background away from the po- 
sition of the source was made and an appropriate amount 
deducted from the measured flux density at the position 



examined. The errors in the corrected flux densities reflect 
the uncertainty in the background level which has been sub- 
tracted. Table ^ shows the measured flux densities of the 
'new' sources from the high-resolution maps and compares 
these measurements with the integrated flux densities from 
the low-resolution maps. Clearly, the flux densities measured 
from the high-resolution maps are systematically lower than 
those measured from the low-resolution maps. This discrep- 
ancy is partly due to the more diffuse parts of the sources 
being 'resolved out' in the 85 mas maps and to the erro- 
neous inclusion of diffuse background in the measurements 
from the 200 mas resolution maps. Figure |E| shows the con- 
tinuum radio spectra of the 20 sources which have not had 
spectra previously published by either W97 or AK98. The 
spectra were fitted by three models for the radio emission 
and 6 examples are shown in Figure ^. The models used can 
be described by the equations 



S(v) = Sou a , 
S(y) = S v a e- r{ "\ 



and 



S(v) = S v 2 (l - 
where, 

r{y) = (v/v T=1 y 



(1) 
(2) 

(3) 

(4) 



These equations represent a simple power-law spectrum, a 
power-law spectrum with free-free absorption by a screen of 
foreground ionised gas and a self-absorbed Bremsstrahlung 
spectrum respectively. From this point they will be referred 
to as models 1, 2 and 3 respectively. In these expressions a is 
the optically thin spectral index and v T= j is the frequency at 
which the free-free optical depth is unity. Of the 20 sources 
which were modelled, 7 were found to be best fitted by model 
1 with an 'inverted' spectral index (a > 0), 1 source was best 
fitted by model 1 with a non-thermal spectral index (a < 
0), 4 sources were best fitted by model 2 and 8 sources by 
model 3. These results are used in the following section to 
make the final source identifications. Note that we interpret 
the spectra of the sources which have inverted spectra with a 
spectral index < 2 (as expected for a simple, optically thick 
Hn region), as being due to an electron density gradient 
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Figure 1. VLA A-configuration and Pie Town image of the central regions of M82 at 15 GHz. The grey-scale range is from to 
1 mjy beam -1 . The dotted lines indicate the fields used to image the MERLIN data and the boxes provide a reference for Figure ^ and 
show the locations within M82 of the spectral index sub-images. Note that this map includes the full spatial frequency range observed 
by the VLA, hence some of this structure has been 'resolved out' in the maps made from the spatial frequency limited data. The boxes 
indicate which sources are shown in Figure ^ due to having significant compact emission at both 5 GHz and 15 GHz and the circles 
indicate the sources which are shown in Figure |l] since they have significant compact emission at 15 GHz only. 



of the type discussed by Olnon (1975 ) and Panagia fc Felli 
(1975|), |since ionised gas with different optical depths which 
are confused will result in an inverted spectral index (for 
the entire optically thick source) which is a function of the 
power- law index of the electron density gradient. 



3.2 Radio source identification 

In order to determine which mechanism is responsible for 
the radio emission from the compact sources in M82, the 
morphology of the sources as resolved in the maps shown 
in Figure g and the spectral information in Figure ^ were 
considered. Therefore, in addition to the 26 which have been 
identified in previous publications, identifications have been 
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Figure 2 - continued 
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Figure 2 - continued 
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made for a further 20 sources. Table ^| lists all of the 46 
sources which have had spectra analysed by W97, AK98 
and in this paper as follows. 

(i) Source identifier. 



(ii) Position of image origin. 



(iii) Optically thin spectral index deduced by Wills et al 



(1997) (W97). 



/) Optically thin spectral index deduced by Allen & 
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Figure 4. VLA A-configuration (incorporating Pie Town) 15 GHz maps of 5 compact sources which remain undetected in the MERLIN 
5 GHz images. 



Table 2. Integrated flux densities of the radio sources for which new spectral information has been presented in this paper. Measurements 
made from the low-resolution (~200 mas) maps and high-resolution (~85 mas) maps are both presented. 



(1) (2) (3) (4) (5) 

Source ID S 5 Sis S 5 Sib 

(85 mas) (85 mas) (200 mas) (200 mas) 



38 


76+534 


0.15 


± 





03 


0.27 


± 





06 


0.26 


± 





05 


0.24 


± 





07 


.39 


29+542 


0.21 


± 





03 


1.3 


± 





1 


1.7 


± 





2 


1.5 


± 





2 


.39 


68+556 


1.1 


± 





1 


3.4 


± 





2 


1.5 


± 





2 


3.2 


± 





2 


40 


62+560 


0.49 


± 





09 


0.22 


± 





05 


1.2 


± 





3 


0.52 


± 





08 


40 


95+588 


0.85 


± 





09 


1.4 


± 





1 


1.4 


± 





4 


2.3 


± 





3 


40 


96+579 










0.62 


± 





07 


0.86 


± 





22 


0.58 


± 





08 


41 


17+562 


0.72 


± 





10 


1.9 


± 





2 


2.0 


± 





3 


2.8 


± 





2 


41 


64+573 


0.36 


± 





06 


1.6 


± 





1 


0.96 


± 





19 


1.8 


± 





2 


42 


08+584 


0.35 


± 





06 


1.6 


± 





1 


1.2 


± 





4 


1.7 


± 





2 


42 


48+584 










0.43 


± 





08 


1.8 


± 





2 


1.0 


± 





2 


42 


56+580 


0.56 


± 





07 


1.5 


± 





1 


2.7 


± 





4 


1.8 


± 





2 


42 


66+564 


0.58 


± 





07 


0.30 


± 





06 


0.62 


± 





25 


0.34 


± 





07 


42 


69+582 


0.54 


± 





06 


1.7 


± 





2 


0.97 


± 





43 


1.8 


± 





2 


42 


82+613 


0.83 


± 





08 


0.41 


± 





07 


1.5 


± 





2 


0.31 


± 





12 


44 


43+618 


0.72 


± 





09 


0.67 


± 





10 


1.1 


± 





3 


0.31 


± 





09 


44 


91+611 


0.71 


± 





07 


0.43 


± 





07 


1.2 


± 





2 


0.39 


± 





11 


44 


93+639 










0.55 


± 





08 


0.26 


± 





09 


0.66 


± 





08 


45 


33+646 










0.59 


± 





09 


0.63 


± 





17 


0.63 


± 





11 


45 


41+637 










0.43 


± 





06 


0.21 


± 





07 


0.48 


± 





10 


46 


17+676 


0.48 


± 





06 


1.0 


± 





1 


1.1 


± 





2 


1.4 


± 





2 
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38.76+53.4 (HII) 





- r 


> 

40.62+56.0 (SNR) 




Frequency (GHz) 



44.91+61.1 (SNR) 



46.17+67.6 (HII) 



Frequency (GHz) 



Figure 5. Example radio continuum spectra of 6 radio sources in M82 which have not had spectra published previously. The lines 
indicate the best fits to the integrated flux densities measured from low-resolution (~0.2") maps. 



Kronberg (1998) - screen absorption model / no absorption 
(AK98 (a)). 



[v] Optic ally thin spectral index deduced by Allen & Kro- 
nberg (1998) - mixed absorption model (AK98 (b)). 



(vi) Two-point spectral index (this paper). 

(vii) Final identification of source as either supernova- 
related (SNR) , Hll region or AGN candidate. 

As can be seen from Figure ^, some of the regions imaged 
contain complex structure, therefore some of the identifica- 
tions listed in Table ^ are more robust than others. For ex- 



ample, the sources at 39.29+54.2, 39.68+55.6, 41.64+57.9, 
42.08+58.4, 42.21+59.2 and 42.56+58.0 as well as the 5 
'new' sources shown in Figure ^ all consist of relatively sim- 
ple structure which is clearly brighter at 15 GHz than 5 GHz. 
Hence, these sources are confidently identified as being the 
most optically thick part of regions of thermal emission. The 
remaining 5 sources which are listed as Hll regions in Ta- 
ble ^ are located in more morphologically complex regions 
which are dominated by thermal processes, although a mix- 
ture of thermal and non-thermal emission is inevitable at 
some level. For example, at 5 GHz, 46.17+67.7 has a simple 
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Table 3. Spect ral indices and identifications of the radio sources which have been deduced by Wills et al. (1997) (W97), 



Allen & 



Kronbei g (1998j ) (AK98) and in this paper (M02). Column (iv) contains the spectral indices deduced by AK98 using a screen absorption 
or no absorption model and column (v) contains the spectral indices deduced by the same authors using a mixed absorption model. 



(i) (ii) (iii) (iv) (v) (vi) (vii) 

Source ID Image Origin a a a a| 5 ID 

(relative to (B1950) 09 h 51 m +69° 54') W97 AK98(a) AK98(b) M02 



38.76+53.4 


38.759 


+53.52 




_ 








_ 




_ 


0.54 ± 0.27 


Hn 


39.10+57.4 


39.106 


+57.33 


-0.6 


± 


0.1 


-0 


38 


± 


0.01 


_ 


-0.53 ± 0.09 


SNR 


39.29+54.2 


39.289 


+54.25 




_ 








_ 




_ 


1.64 ± 0.16 


Hn 


39.40+56.1 


39.391 


+56.18 


-0.50 


± 


0.05 


-0 


21 


± 


0.01 


_ 


-1.04 ± 0.22 


SNR 


39.64+53.4 






-0.20 


± 


0.05 


-0 


71 


± 


0.01 


_ 


_ 


SNR 


39.68+55.6 


39.668 


+55.59 




_ 








_ 




_ 


1.03 ± 0.09 


Hn 


39.77+56.9 






-0.50 


± 


0.06 


-0 


19 


± 


0.01 


_ 


_ 


SNR 


40.32+55.1 


40.317 


+55.20 


-0.50 


± 


0.06 


-0 


55 


± 


0.01 


_ 


-0.23 ± 0.21 


SNR 


40.62+56.0 


40.594 


+56.10 
















_ 


-0.72 ± 0.25 


SNR 


40.66+55.2 


40.676 


+55.11 


-0.70 


± 


0.05 


-0 


52 


± 


0.01 


-0.55 ± 0.01 


-0.54 ± 0.08 


SNR 


40.95+58.8 


40.938 


+58.87 




_ 








_ 




_ 


0.44 ± 0.13 


Hn 


40.96+57.9 


40.955 


+57.85 




_ 








_ 




_ 


> 1.2 


Hn 


41.17+56.2 


41.176 


+56.16 




_ 














0.87 ± 0.14 


Hn 


41.29+59.7 


41.302 


+59.64 




_ 




-0 


54 


± 


0.01 


-0.56 ± 0.01 


-0.47 ± 0.09 


SNR 


41.64+57.9 


41.637 


+57.94 




_ 








_ 




_ 


1.32 ± 0.18 


Hn 


41.95+57.5 


41.951 


+57.49 


-0.80 


± 


0.05 


-0 


72 


± 


0.01 


-0.74 ± 0.01 


_ 


SNR? 


42.08+58.4 


42.079 


+58.39 




_ 












_ 


1.32 ± 0.17 


Hn 


42.21+59.2 


42.210 


+59.04 




_ 




-0 


10 


± 


0.01 




1.16 ± 0.13 


Hn 


42.48+58.4 


42.481 


+58.36 




_ 








_ 




_ 


> 1.2 


Hn 


42.53+61.9 






-0.80 


± 


0.05 


-1 


84 


± 


0.01 


-1.99 ± 0.01 


_ 


SNR 


42.56+58.0 


42.557 


+58.05 
















_ 


0.88 ± 0.14 


Hn 


42.66+56.4 


42.673 


+56.40 


















-0.58 ± 0.21 


SNR 


42.67+55.6 


42.662 


+55.54 


-0.30 


± 


0.05 


-0 


61 


± 


0.01 


-0.63 ± 0.01 


-1.3 ± 0.2 


SNR 


42.69+58.2 


42.694 


+58.24 




_ 








_ 




_ 


1.04 ± 0.13 


Hn 


42.82+61.3 


42.810 


+61.30 




_ 








_ 




_ 


-0.63 ± 0.17 


SNR 


43.18+58.3 


43.186 


+58.35 


-0.80 


± 


0.05 


-0 


67 


± 


0.01 


-0.72 ± 0.01 


-0.44 ± 0.08 


SNR 


43.31+59.2 


43.305 


+59.20 


-0.60 


± 


0.06 


-0 


61 


± 


0.01 


-0.85 ± 0.01 


-0.65 ± 0.07 


SNR 


44.01+59.6 


44.008 


+59.59 


0.20 


± 


0.05 


-0 


51 


± 


0.01 


-0.56 ± 0.01 


-0.38 ± 0.06 


SNR 


44.29+59.3 


44.276 


+59.30 








-0 


56 


± 


0.01 




-0.72 ± 0.12 


SNR 


44.43+61.8 


44.419 


+61.72 


















0.07 ± 0.17 


SNR 


44.52+58.1 


44.509 


+58.21 


-0.60 


± 


0.05 


-0 


61 


± 


0.01 


-0.64 ± 0.01 


-0.15 ± 0.17 


SNR 


44.91+61.1 


44.899 


+61.19 


















-0.45 ± 0.18 


SNR 


44.93+63.9 


44.934 


+63.91 


















> 1.1 


Hn 


45.17+61.2 


45.173 


+61.25 


-0.80 


± 


0.05 


-0 


68 


± 


0.01 


-0.69 ± 0.01 


-0.52 ± 0.07 


SNR 


45.26+65.3 


45.254 


+65.18 








-0 


62 


± 


0.01 




-1.02 ± 0.20 


SNR 


45.33+64.6 


45.330 


+64.63 


















> 1.1 


Hn 


45.41+63.7 


45.409 


+63.73 


















> 1.0 


Hn 


45.44+67.3 


45.423 


+67.43 


-0.60 


± 


0.05 


-0 


57 


± 


0.01 




-0.83 ± 0.17 


SNR 


45.52+64.7 






-0.20 


± 


0.09 


-0 


15 


± 


0.01 






SNR 


45.79+65.2 


45.758 


+65.32 


-0.60 


± 


0.05 


-0 


23 


± 


0.01 




-0.55 ± 0.13 


SNR 


45.91+63.8 


45.898 


+63.88 


-0.20 


± 


0.05 


-0 


53 


± 


0.01 


-0.54 ± 0.01 


-0.38 ± 0.11 


SNR 


46.17+67.6 


46.175 


+67.72 


















0.66 ± 0.16 


Hn 


46.52+63.8 


46.521 


+63.93 


-0.40 


± 


0.05 


-0 


73 


± 


0.01 


-1.50 ± 0.01 


-0.35 ± 0.11 


SNR 


46.56+73.8 






-0.60 


± 


0.05 


-0 


78 


± 


0.01 


-1.04 ± 0.01 




SNR 


46.75+67.0 






-0.80 


± 


0.05 


-0 


57 


± 


0.01 


-0.90 ± 0.01 




SNR 


47.37+68.0 






-0.60 


± 


0.07 


-0 


57 


± 


0.01 






SNR 



single-component structure which becomes more complex 
at 15 GHz. A mixture of thermal and non-thermal emis- 
sion may also be responsible for the structures of sources 
which have been identified as supernova remnants, exam- 
ples of which are 44.43+61.8 and 44.52+58.2. 

Of the 26 sources identified previously, only one was 
deemed to be a possible Hn region. This number has now 
increased to 16 by the addition of 15 sources examined in 
this paper. In addition, a further 5 supernova remnants have 
been identified. Therefore, it is suggested that the previous 



studies were biased towards the identification of supernova 
remnants, since only the most luminous sources were in- 
vestigated. Also, the measurements of the source spectra 
were also biased towards longer wavelengths, away from the 
regime where thermal sources are brightest - especially if 
they are optically thick. The observations made with the 
VLA and Pie Town combination have provided the highest 
resolution images of the central regions of M82 at 15 GHz 
and clearly allowed an extension of previous work into the 
thermal regime. The peak brightness temperature of a radio 
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■ SUPERNOVA REMNANTS (established) 
□ SUPERNOVA REMNANTS (newly identified) 

• HII REGIONS (established) 
O HII REGIONS (newly identified) 




Figure 6. Two-point spectral index versus brightness tempera- 
ture plot for 46 compact radio sources in M82. 



source can also be used as a constraint when identifying Hn 
regions. The brightness temperature, Tt, is related to the 
electron temperature, T e by 



T b = T e (l - e~ T ), 



(5) 



where r is the optical depth. Hence in the optically-thick 
regime (r 3> 1) Tt ~ T e . Since ionised gas becomes opaque 
towards lower frequencies it was decided to measure the 
brightness temperatures from the lowest frequenc y map 
availab le, which w as the L-band (1420 MHz) map of Pedlar 



& Mux ow (1995 ) and Figure ^ shows the measured bright- 



ness temperatures plotted against the source spectral in- 
dices. It can be seen that the brightness temperatures of the 
Hn regions all lie below 10 4 K, as indicated by the horizontal 
line. It should be noted that the measured brightness tem- 
perature of a thermal source is a lower limit for the electron 
temperature because as an Hn region becomes optically thin 
the brightness temperature becomes less than the electron 
temperature. It can also be seen that the brightness tem- 
peratures derived for the supernova remnants are mostly in 
excess of 10 4 K, although since the previous studies of the 
spectra by W97 and AK98 concentrated only on the bright- 
est sources this can be mainly attributed to a selection effect. 



3.3 The radio source distribution in M82 

Since the number of identified sources in M82 is now 46, 
a rudimentary statistical analysis of the distributions may 
be attempted. Comparing the distribution of Hn regions 
with identified supernova remnants, it can be said that 
the two distributions are different with 80% confidence. 
This confidence level was der ived using the two-dimensiona l 
Kol mogorov-Smirn ov test of Fasano & Franceschini (1987) 
and Peacock (1983). Figure ^ shows the one-dimensional dis- 
tributions of the right ascension coordinates of the compact 
radio sources in M82. The top panel shows the distribution 
of Hn regions and the bottom panel shows the distribution of 
supernova remnants, including an add itional 14 un identified 
sources. These sources were noted by Wills (1998) but were 
undetected at 408 MHz and were not included in the spectral 
analysis in W97. However, these sources are almost certainly 




SUPERNOVA REMNANTS 



Right Ascension (relative to (B1950) 09 h 51" | 



Figure 7. The one-dimensional distributions of the compact ra- 
dio sources in M82 projected onto the right ascension axis. The 
top panel shows the distribution of Hn regions and the bottom 
panel shows the distribution of supernova remnants. 



supernova remnants, since they are generally present on low- 
frequency (~1.6 GHz) maps, but not at higher frequencies - 
implying a non-thermal spectrum. Due to their low bright- 
nesses a spectral analysis has yet to be made for them. 
In addition, by taking into account the transient source at 
41.50+59.7, the total number of compact sources whose ori- 
gin has been determined is now 61. The propagatio n of star- 
forming regi ons w ithin M82 has been d iscussed by satyapal 



ct al. (1997) and de Grijs et al. (2000). On the basis of a 



changing CO index across M82 and other indicators of an 
old stellar population in the nucleus, Satyapal et al. (1997) 
suggested that the st arburst was propagat ing outwards from 
the centre. However, de Grijs et al. (200C) reported the dis- 
covery of 'fossil' supernova remnants away from the nucleus 
(~ 500 pc away) and suggested that the star-formation was 
propagating inwards. In order to test these hypotheses with 
the radio data, the sizes of the Hn region and SNR distribu- 
tions shown in Figure [?] was examined, since the Hn regions 
trace the hot, young stars and are more indicative of re- 
cent star-formation than the supernova remnants which are 
the result of star- formation ~10 7 years previously. Figure ^ 
shows how the positions of the 16 Hn regions correlate with 
the large scale distribution o f ionised gas in M82 as sho wn 
by the 12.8 /im [Nell] map of Achtermann fc Lacy (1991: ). It 
can be seen that the [Nell] distribution has a width of ~20" 
at half-intensity which compares with a width of ~35" for 
the supernova remnant distribution in Figure [j] which may 
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Table 4. Properties of the Hn regions in M82. 



RIGHT ASCENSION (B1950) 



Figure 8. The ionised gas in M82. The crosses mark the posi 
tions of the compact Hll regions identified i n this paper and are 



over l aved on the [Nell] line distribution of Achtermann & Lacy 
(1993). [The contour levels are at integer multiples of 10% of the 
peak intensity. 



Source 


Size 


Emission Measure 


Equivalent 


ID 


(PC) 


(cm -6 pc) 


number of 








(xlO 7 x (Te/lO 4 ^) 1 - 35 ) 


05 stars 


38.764 


-53.4 


1 9 


1 2 


3 


39.294 


-54.2 


1 8 


2 5 


5 


39.684-55.6 


5 6 


21 


520 


40.954 


-58.8 


7.2 


13 


530 


40.964-57.9 


2.3 


2.0 


8 


41.174 


-56.2 


7.9 


10 


490 


41.644 


-57.9 


5.6 


>9.0 


221 


42.084 


-58.4 


4.3 


5.9 


86 


42.214 


-59.2 


1.9 


5.0 


94 


42.484 


-58.4 


0.8 


3.4 


2 


42.564 


-58.0 


4.3 


2.3 


33 


42.694 


-58.2 


1.9 


>9.0 


25 


44.934 


-63.9 


3.7 


>9.0 


97 


45.334 


-64.6 


2.0 


1.3 


4 


45.414 


-63.7 


1.3 


>9.0 


12 


46.174-67.6 


3.6 


6.9 


70 



seem to fav our the inwardly prop agating; star-formation hy- 
pothesis of de Grijs et al. (200C). However, small-number 



statistics in the compact source distribution do not allow a 
robust distinction to be drawn between the two possibilities. 

Again referring to Figure it can be seen that a num- 
ber of the thermal sources appear to be associated with the 
two peaks in the [Nell] emission in the west, although a 
number of the sources seem to be in regions of low [Nell] 
intensity. However, another two-dimensional Kolmogorov- 
Smirnov test comparing the compact Hll region distribution 
with the [Nell] line distribution determined that the distri- 
butions are similar with a significance level of 90%. 



3.4 Compact radio sources of non-thermal origin 

For the 30 non-thermal sources, it can be seen from Table ^ 
that the measured two-point spectral indices are often signif- 
icantly different from those derived for the models of W97 
and AK98. As shown in Figure ^, these sources are often 
highly resolved and more morphologically complex than the 
compact Hll regions identified previously. It would appear 
that the two-point spectral indices are more consistent with 
the models for the brightest, most compact sources and less 
consistent for the weaker, complex sources. This discrepancy 
may be partly due to the 'resolving out' of larger scale emis- 
sion. In any case, deciding on what size scale some of the 
sources 'start' is a very subjective matter. 



3.5 The compact radio sources of thermal origin 

In the previous studies of the spectra of the most compact 
radio sources in M82, only two sources have been identi- 
fied as being of possible thermal origin, 40.624-56.0 (W97) 
and 42.214-59.2 (AK98). It should be noted at this stage 
that we identify 40.624-56.0 as a supernova remnant due 
to the extension of the spectrum of this source to higher 
frequencies than those presented in W97. The spectrum of 



this source is shown in Figure ^ and if the 15 GHz measure- 
ment is excluded then the spectrum agrees with that in W97 
and appears to be flat. The higher frequency measurement, 
however, implies a non-thermal orig in for the radio emis- 
sion from this source. In Section 3.2, 16 sources (including 
42.214-59.2) were identified for which positive or flat spectral 
indices have been derived which are attributable to free-free 
emission from ionised gas. Table ^ summarises the inferred 
properties of these Hll regions. For the resolved sources, the 
sizes were measured by taking a slice through the emission 
and determining the largest angular size. For the more com- 
pact sources, a two-dimensional Gaussian was fitted to the 
source using the AIPS task jmfit. The emission measures 
were derived from the fits to the spectra of the sources as 
described in Section [Tl]. For the cases where the sources 
did not display a turnover in their spectra at the transi- 
tion between the optically thick and optically thin regimes, 
a lower limit for the emission measures was taken to be 9 x 
10 7 cm -6 pc. This limit corresponds to a turnover frequency 
somewhere above 5 GHz. By calculating the Lyman contin- 
uum flux required to ionise the Hll regions, given their sizes 
and emission measures stated in Table ^ it is possible to 
derive an equivalent number of 05 stars, assuming that the 
Lyman co ntinuum out put of an 05 star is ~5.1 x 10 49 pho- 
tons s _1 (Bpitzer 1975). The number of 05 stars required to 
ionise these regions varies from ~ a few, to over 500. There- 
fore, these regions may be similar to those fo und in other 
galaxies such as NGC 2146 (jTarchi et al. 2000| ), Henize 2-10 
( |Kobulnicky fc Johnson 1999[ ) and NGC 5253 (Turner, Beck 
& Ho 2000) , which were deemed to be due to the precursors 
of 'super-star clusters' or 'proto-globular clusters'. No opti- 
cal counterparts were found for the radio sources in these 
galaxies, however, since these galaxies are not face-on, the 
optical emission is likely to be obscured by intervening dust 
and/or the parent molecular cloud. 

With the relatively high resolution provided by the VLA 
and MERLIN (85 mas = 1.3 pc) it is now possible to resolve 
the compact thermal sources. Variations in the two-point 
spectral indices can be seen across several of the sources 
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For the simple model of an Hn region in which S oc v in 
the optically-thick regime and the optically-thin 

regime this would imply that these sources are becoming op - 



GAUSSIAN MODEL 

POWER LAW MODEL 

CONSTANT DENSITY SPHERE MODEL 



Figure 9. Two-point spectral index (between 5 and 15 GHz) vari- 
ations across the Hn region at 42.21+59.2. The lines indicate the 
best fits to the profile for 3 different electron density distributions. 



Table 5. The number of identified compact Hn regions and young 
supernova remnants in M82. 





Before 


This work 


TOTAL 




this work 






Supernova Remnants 


25 


5 


30 


Hn Regions 


1 


15 


16 



and Figure ^| shows the spectral index variation across the 
brightest Hn region at 42.21+59.2. Three models for the 
electron density distribution were fitted to these data - a 
power-law, Gaussian distribution and constant electron den- 
sity sphere. As can be seen from Figure || it is currently not 
possible to choose between these models. However, the in- 
creased spectral index towards the centre of the region may 
be attributed to increasing optical thickness due to enhanced 
emission measure. 



4 CONCLUSIONS 

The VLA in its A-configuration connected with the Pie 
Town VLBA antenna via a real-time fibre optic link has 
been used to image M82 at 15 GHz. These observations have 
provided the highest resolution images yet made at this fre- 
quency and allowed a direct comparison with 5 GHz MER- 
LIN images at a similar resolution. The conclusions may be 
summarised as follows. 

(i) In addition to the 26 compact radio sources which have 
been previously identified in the nucleus of M82, a further 
20 sources have been identified. The final identifications are 
summarised in Table [| (n.b. the controversial sources at 
41.95+57.5 and 44.01+59.5 are included under 'supernova 
remnants'). By comparison with the 12.8 /im [Nell] map of 
Achtermann fc Lacy (1995 ), it was noted that the locations 
of the compact Hn regions broadly coincided with the peaks 
in the large-scale ionised gas distribution. 

(ii) Of the 15 new Hn regions identified, all have steep 



tic al ly thick between 5 and 15 GHz. However, Olnon (197! 
and Panagia & Felli (1975) showed that simple models of 
ionised nebulae with significant electron density gradients 
could produce inverted spectra, the slope of which depends 
on the structure of the region. This may be the case for a 
number of the sources, but further high-resolution imaging 
at different frequencies is required to confirm or deny this 
possibility. 

(iii) The low-resolution spectra of the regions in which the 
compact sources have been imaged show inverted spectra 
for about half of the sources and flattening spectra at the 
highest frequencies for the remainder. It is inferred that at 
the highest resolution only the most compact and optically- 
thick components of the ionised gas are being detected. 

(iv) Of the 46 identified sources in M82, ~ 35% are Hn 
regions and the remaining ~ 65% probably have supernova 
origins. Therefore, it would appear that M82 is not lacking 
in compact H n regions in compa rison to other galaxies such 
as NGC 2146 (Tarchi et al. 200C) as much as was previously 



thought, although it still appears that M82 is in a more 
advanced starburst stage than NGC 2146 due to the greater 
number of supernova remnants. 

(v) Five additional sources were identified as supernova 
remnants since they have steep, negative spectral indices 
at the higher frequencies. At the lower frequencies they are 
often experiencing a turnover in their spectra by 1.4 GHz 
requiring emission measures for the foreground ionised gas 
of ~ 10 7 cm -6 pc. These supernova remnants are obscured 
at the lower radio frequencies where previous studies have 
concentrated and may represent those supernova remnants 
embedded deeper within M82 and hence behind a greater 
amount of ionised gas. 

(vi) An additional 14 unidentified sources were added to 
the analysis and after taking into account the transient 
source noted by Kronberg, Biermann & Schwab (1985), the 
total number of compact sources now stands at 61. The 14 
unidentified sources are most likely the older supernova rem- 
nants in the sample and hence are not bright enough to have 
been included in the essentially flux-limited analyses. How- 
ever, assuming that they are all SNR, a comparison between 
the Hn region distribution and SNR distribution shows sig- 
nificant differences, although any difference in the overall 
extent of the distributions is not significant. Hence, the ra- 
dio data cannot be used to either con firm or deny the prop - 
agat ing star-formation hy potheses of Batyapal et al. (1997) 
and|de Grijs et al. (200of). 



'inverted' spectra at the highest resolution (a 5 > +0.4). 
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